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Abstract: A complete parton level analysis of £'^i~ + four jets {i = e, fi) and Sii' + 
two jets production at the LHC is presented, including all processes at order 0{a%j^j), 
0(a|J^^a|) and 0{a%j^,jag) when appropriate. The infinite Higgs mass scenario, which 
is considered as a benchmark for strong scattering theories and is the limiting case for 
composite Higgs models, and one example of a model incorporating a Strongly Interacting 
Light Higgs are confronted with the Standard Model light Higgs predictions. This analysis 
is combined with the results in the iu + four jets channel presented in a previous paper, 
in order to determine whether a composite Higgs signal can be detected as an excess of 
events in boson-boson scattering. 
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1. Introduction 

The Standard Model (SM) describes Electroweak Symmetry Breaking (EWSB) in the most 
simple and economical fashion through a single complex Higgs doublet, with only a neutral 
scalar field in the spectrum. The fit of EW precision data is in agreement with the SM 
predictions to an unprecedented accuracy and gives an upper limit on the Higgs mass 
of about 182 GeV [1]. Direct searches at LEP2 imply niH > 114 GeV [2] while, more 
recently DO and CDF at the Tevatron have excluded at 95%CL a SM Higgs in the range 
160 GeV < TTiH < 170 GeV [3]. The LHC will have the task to reveal whether this minimal 
realization of EWSB takes place in Nature or a more complex structure is present^. 

If the Higgs field is not found, scattering processes between longitudinally polarized 
vector bosons will play a prominent role because, without a Higgs, the corresponding 
amplitudes grow with energy and violate perturbative unitarity at about one TeV [9], 
requiring new physics in the energy range accessible to the LHC. 

Many alternative mechanisms of EWSB have been explored. We will not try to sum- 
marize the different models and simply refer to the literature. For our purpose we will 
only remark that it is conceivable and widely discussed [10, 11, 12, 13, 14, 15, 16] that 
composite states are responsible for EWSB as nicely recently reviewed in Ref. [17]. These 

^Detailed reviews and extensive bibliographies can be found in Refs. [4, 5, 6, 7, 8] 
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theories are characterized by the presence of new states which could be produced at the 
LHC, if hght enough. 

In view of the large number of different proposals it is useful to determine the model 
independent features of this class of theories. There has been recent progress in this area 
[18, 19, 20], using the effective theory language [21]. In Ref. [19] it has been pointed out 
that, if EWSB is triggered by a light composite Higgs which is a pseudo-Goldstone boson 
related to some large scale strongly interacting dynamics, the growth with energy of the 
vector boson scattering amplitudes typical of Higgsless models might not be completely 
canceled by Higgs exchange diagrams but only slowed down because of the modified cou- 
plings between the Higgs field and the vector bosons with respect to the SM ones. Such a 
Higgs has been called Strongly Interacting Light Higgs (SILH). If the mass of the Higgs bo- 
son becomes larger than the typical energy scale at which boson-boson scattering is probed, 
the contribution of the Higgs exchange diagrams decreases and completely vanishes, in the 
Unitary gauge we employ throughout our calculation, in the limit of an infinite mass which 
we will refer to as the Higgsless case. As a consequence, the scattering cross section for an 
infinitely massive Higgs in the SM represents, at large energies, an upper limit for VV scat- 
tering processes in SILH models, and can be taken as a benchmark for the observability of 
signals of strong scattering and of Higgs compositeness in boson-boson reactions. On the 
other hand the Higgsless case is also representative of models in which heavy resonances 
which unitarize boson-boson scattering are present but cannot be directly detected at the 
LHC. 

Scattering processes among vector bosons have been scrutinized since a long time 
[22, 23]. In Ref. [24, 25] an analysis of + four jets and + four jets produc- 

tion at the LHC has been presented, with the limitation of taking into account only purely 
electroweak processes. Preliminary results concerning the inclusion of the 0{a%,^a1) back- 
ground, which include VV + 2j and top-antitop production have appeared in Ref. [26]. A 
preliminary analysis in the Equivalent Vector Boson Approximation of the observability of 
partial unitarization of longitudinal vector boson scattering in SILH models at the LHC 
can be found in Ref. [27]. In the last few years QCD corrections to boson-boson production 
via vector boson fusion [28] at the LHC have been computed and turn out to be below 
10%. Recently, VBFNLO [29] a Monte Carlo program for vector boson fusion, double and 
triple vector boson production at NLO QCD accuracy, limited to the leptonic decays of 
vector bosons, has been released. 

In Ref. [30] a complete parton level analysis of iv + four jets production at the LHC, 
including all processes at order 0{a%j^), 0{a%j^ja^^) and 0{a%j^ag) has been presented, 
comparing a typical SM light Higgs scenario with the Higgsless case. It was noted that 
the 0{a%j^ag) W + A] background is so large to make the usual approach of comparing 
the number of events in the two scenarios at large invariant masses rather ineffective. An 
alternative strategy was suggested, namely to focus on the invariant mass distribution of 
the two central jets which is characterized by a peak corresponding to the decays of vector 
bosons which is present in the vector-vector scattering signal and by an essentially flat 
background produced by C'(a|J^,^Q|) W +Aj processes. The latter can be measured from the 
sidebands drastically decreasing the theoretical uncertainties. In Ref. [30] the probability 
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of observing a signal of new physics in the W + 4j channel was estimated constructing 
the probability distribution of experimental results both in the SM and in the Higgsless 
case, taking into account the statistical uncertainties for both signals and background and 
the theoretical uncertainty for the signal alone, assuming that the background can be 
extrapolated from the measured sidebands. The probability that the benchmark Higgsless 
scenario would result in an experimental outcome outside the SM 95% exclusion limit, 
assuming a 200 fb~^ luminosity for both the muon and electron decay channel, turned out 
to be about 97%. 

In this paper we complete our study examining at parton level the processes pp — ^ 
£'^i~ + 4j and pp Siu + 2j, including all irreducible backgrounds contributing to these 
six parton final states. These are the only vector-vector scattering reactions, together 
with pp ^ iu + 4j mentioned above, in which the presence of at least one pair of leptons 
in the final state ensures a sufficiently clean determination of the invariant mass of the 
outgoing vector bosons, up to the ambiguity related to the momentum reconstruction of 
the unobserved neutrino. The latter can be however be estimated by imposing the condition 
that the missing momentum and the charged lepton reconstruct the W mass. It is here 
implied that the pp ^ 41 + 2j channel, while allowing a more precise reconstruction of the 
VV mass, has a negligible rate at the LHC at large invariant masses of the four- lepton 
system due to the small branching ratio of the leptonic Z decays. The invariant mass 
Mvv of the two outgoing bosons is the analogue of the center-of-mass energy for on-shell 
vector boson scattering: differences between the Higgsless scenario and the SM predictions 
increase, for the processes we analyze, with the boson pair invariant mass. 

In this paper we consider three scenarios: a light Higgs SM framework with Mh = 200 
GeV, one instance of the SILH models which we will describe shortly and an infinite mass 
Higgs scenario. We combine the results for pp i^i^ + 4j and pp 2>lv + 2j with those 
obtained in Ref. [30] for the pp ^ iv + 4j channel in order to obtain a global estimate, 
including all three channels, of the probability of distinguishing the considered Beyond 
Standard Model (BSM) cases from the SM scenario. 

2. The Strongly Interacting Light Higgs scenario 

The effective field theory approach [21] is a powerful tool for describing the low energy 
dynamics of systems with broken symmetries. It provides a systematic expansion of the 
full unknown Lagrangian in terms of the fields which are relevant at scales much lower than 

the symmetry breaking scale. In Ref. [19] this general framework has been employed to 
describe the main features of the models of EWSB in which the Higgs field can be identified 
with a pseudo-Goldstone boson of some broken strong interaction at high energy. Models 
which fall into this class are for instance the Holographic Higgs [15], the Little Higgs of 
Ref. [16] and the Littlest Higgs [12]. 

These theories generally predict new resonances and their search may well be the best 
way to detect New Physics effects at the LHC. However they also lead to modifications 
of boson-boson scattering which might provide further evidence and possibly be the only 
accessible effects if the new resonances are heavy. 
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The leading low energy effects are described by two parameters (one responsible for a 
universal modification of all Higgs couplings, and the other one for a universal modification 
of Higgs couplings to fermions) characterized by the ratio / p = ^, where v is the Higgs 
vacuum expectation value and / is the cr-model scale. The natural range of the ^ parameter 
is between ^ = and C = 1 which correspond respectively to the limiting cases of the 
Standard Model and of technicolor theories. Because of the modified Higgs couplings, 
longitudinal gauge-boson scattering amplitudes violate unitarity at high energy, even in 
the presence of a light Higgs [19]. The i?^-growing amplitude is typically smaller than in 
the Higgsless case and the violation is moved to a larger energy regime. 

Therefore, even if a light Higgs is discovered, boson-boson scattering is a crucial pro- 
cess to study and can give us useful information on the nature of the Higgs boson. It 
is worth pointing out that, in this framework, since the Higgs can be viewed as an ap- 
proximate fourth Goldstone boson, its properties are related to those of the exact (eaten) 
Goldstone bosons. Strong gauge-boson scattering will be accompanied by strong Higgs 
pair production [19]. 

The effective Lagrangian approach of Ref. [19] is valid for small values of ^, while 
larger values demand a more detailed description of the particular model at hand. Such a 
Lagrangian leads to a modification of the Higgs couplings by a factor 1/-^1 + c//^, which 
can be reabsorbed in a Higgs propagator modification by a factor 1/(1 + chC) boson 
boson scattering studies, ch is a pure number of order unity [12, 15, 16, 19]. For the 
present study we have selected the value ch(, = 1 which we intend as a possible upper limit 
for the model independent lagrangian description of Ref. [19] . 

3. Outline of the analysis 

The observation of strong boson-boson scattering as an excess of events compared to the 
SM prediction requires, as an essential condition, that a signal of VY scattering is extracted 
from the background. At the same time the selection strategy must be capable to max;imize 
the differences between the light Higgs and the BSM cases. 

Three(two) perturbative orders contribute to the background for + 4j{3ii' + 2j) 
final states. At 0{a%^) there are a large number of diagrams which cannot be interpreted 
as boson-boson scattering and which cannot be separated in any sensible way from the 
scattering type diagrams due to large cancellations between the two sets [25]. For illus- 
trative purpose a number of diagrams are shown in Figs. 1-6. Fig. 1 illustrates the class 
of diagrams which include boson-boson scattering subdiagrams while Fig. 2 presents some 
representative 0(a|^^) diagrams in which two vector bosons are produced but no fusion 
takes place. Fig. 3 and Fig. 4 show diagrams which describe Higgs production and three 
vector boson production. 

At 0{a%i^a^) we have to deal with the production of two electroweak bosons plus two 
jets without any scattering contribution. Some representative diagrams are shown in Fig. 5 
which also includes diagrams describing fl production which contribute only to the iv + 4j 
channel. For the + 4j final state there are contributions at O^a^^j^jag) in which only 
one electroweak boson is effectively produced, while the additional jets, which, taken in 
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pairs, do not peak at any particular mass, populate the full available phase space with a 
production rate which is much larger than the signal one. A few diagrams in this set are 
presented in Fig. 6. 

The first step is concerned with the identification of a suitable kinematic signature 
which allows to capture the essence of VV scattering. The selection of events widely 
separated in pseudorapidity is a well established technique for enhancing the scattering 
contributions at the LHC [22, 23]. Looking at the topology of the diagrams embedding the 
gauge boson scattering as a subprocess, one concludes that it is appropriate to associate the 
two most forward/backward jets to the tag quarks which radiate the bosons which initiate 
VV scattering. As shown in Ref. [30] a powerful tool to increase the separation between the 
SM predictions and those of the Higgsless scenario is provided, at large invariant masses, 
by the request that the vector bosons and their decay products are in the central part 
of the detector since the vector bosons in the Higgsless case have smaller rapidities and 
larger momenta than in the presence of a light Higgs. The main purpose of this kinematic 
selection is to isolate a sample of genuine VV+2j events while suppressing the contribution 
of irreducible backgrounds such as three boson production or top quark production. 

Having isolated a sample of candidate scattering events, one needs to define an ob- 
servable quantity which is as susceptible as possible to the details of the mechanism of 
EWSB in order to maximize the sensitivity to effects of alternative models such as strong 
scattering. For the i^i^ + 4j channel this task is straightforward only apparently. As 
already mentioned, the QCD background is expected to provide a significant number of 
fake VV + 2j events as a consequence of the large cross section. The classical approach is 
to focus on the invariant mass distribution of the final state boson pair. The huge QCD 
background, with its large scale uncertainty, makes this procedure rather dubious for the 
+ 4j case. Detecting signals of EWSB in the vector pair mass distribution remains a 
formidable challenge. The analysis of Ref. [30] shows that, even focusing the attention on 
the high invariant vector pair mass region and requiring the mass of the two central jets to 
be close to the mass of a vector boson, the expected signal over background ratio is of the 
order of 1/10. This way of measuring the signal, typical of a counting experiment, finds a 
substantial obstacle in the spread of the interesting excess of events over a wide region in 
Mw 

A possible way out for this problem is to look instead at the invariant mass of the two 
central jets (Mj^j^) for events with large vector pair mass. Provided a convenient set of 
kinematic cuts has been applied, the 0{a%j^) + 0{a%j^a^) cross section is dominated by the 
peaks corresponding to W and Z decays to quarks, while the 0{a'^j^jaj,) Z+Aj contribution 
is non-resonant in this respect. When restricting to the window between 70 and 100 GeV, 
which covers completely the W and Z resonances, we find that the Z + 4j distribution is 
essentially fiat and therefore can reliably be measured from the sidebands of the physical 
region of interest. This procedure has several advantages. On one side, it eliminates the 
theoretical uncertainty associated with the scale dependence of the 0{a%j^a^) contribution. 
On the other side, it allows to subtract the dominant contribution to the irreducible QCD 
background, enhancing the visibility of genuine EWSB effects. The signal has a very clear 
signature, a bump in the two central jet mass distribution, which is much easier to hunt 



-5- 



q q q q q q 




q q q q q q 



Figure 1: Vector boson fusion processes. 




Figure 2: Two vector boson production witliout fusion. 



for experimentally than a diffuse excess in production cross section. 

Once the non-resonant background has been subtracted, one is left with a peak whose 
size is strictly related to the regime of the EWSB dynamics: a strongly-coupled scenario 
would result in a more prominent peak than a weakly-coupled one. This feature suggests 
to take the integral of the peak as the discriminator among different models. At a given 
collider luminosity, the number of expected events can be derived. With a slight abuse of 
language, we call this number the VBS signal. It is by analyzing the probability density 
function (p.d.f.) associated with this discriminator that we can determine, in the last step, 
the confidence level for a given experimental result to be or not to be SM-like, in the same 
spirit of the statistical procedure adopted for the search of the Standard Model Higgs boson 
at LEP [31]. 

For the 3£i/ + 2j channel the physical picture is straightforward, the invariant mass of 
the vector boson pair can be directly measured from the momenta of the charged leptons 
and the reconstructed neutrino momentum. Since the background is much smaller than 
in the + Aj case we use the total expected number of events as discriminator and 

estimate the probability that the experimental results in the BSM models are incompatible 
with the SM. 

4. Calculation 

As discussed in Sect. 3, three perturbative orders contribute to +Aj at the LHC, while 
only two perturbative orders contribute to Zlu + 2j. The 0[a%,^j) and ©(Q^j^^al) samples 
have been generated with PHANTOM, a dedicated tree level Monte Carlo generator which is 
documented in Ref. [32] while additional material can be found in Refs. [33, 34, 35]. The 
0{a\j^ja'^) sample has been produced with MADEVENT [36]. Both programs generate events 
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Figure 3: Higgs boson production via vector boson fusion and Higgsstrahlung. 




Figure 4: Three vector boson production. 




Figure 5: Examples of contributions to the QCD irreducible background: VV + 2j 



in the Les Houches Accord File Format [37]. In all samples full 2 — > 6 matrix elements, 
without any production times decay approximation, have been used. For the LHC we have 
assumed the design energy of 14 TeV. For each perturbative order we have generated a 
sample of five hundred thousand unweighted events. 

For the Standard Model parameters we use the input values: 

Mw = 80.40, Mz = 91.187 GeV, 

Gf, = 1.16639 10-^ GeV-2, as(Mz) = 0.118 

Mt = 175.0 GeV, Mb = 4.8 GeV. (4.1) 

The masses of all other partons have been set to zero. We adopt the standard G^- 
scheme to compute the remaining parameters. 

All samples have been generated using CTEQ5L [38] parton distribution functions. 
For the 0{a%j^f) and 0{a%j^ja1) samples, generated with PHANTOM, the QCD scale has been 
taken as: 

1 ^ 

Q^ = M^ + -Y,pli (4.2) 
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(a) (b) (c) (d) 

Figure 6: 

Representative Feynman diagrams for the 0{a1,j^,ja'^), V + 4j production processes at the 

LHC. 

while for the C'(a|j,,^a|) sample the scale has been set to = M^. This difference in 
the scales conservatively leads to a definite relative enhancement of the Z + 4j background. 
Tests in comparable reactions at 0(a|j^^ag) have shown an increase of about a factor of 
1.5 for the processes computed at = M| with respect to the same processes computed 
with the larger scale Eq.(4.2). 

We work at parton level with no showering and hadronization. The two jets with 
the largest and smallest rapidity are identified as forward and backward tag jet 
respectively. The two intermediate jets called central (jc) and are considered as candidate 
vector boson decay products. 

The neutrino momentum is reconstructed according to the usual prescription, requiring 
the invariant mass of the lu pair to be equal to the W boson nominal mass, 

(/+P-)2 = m2., (4.3) 

in order to determine the longitudinal component of the neutrino momentum. 

For very large Higgs masses, all Born diagrams with Higgs propagators become com- 
pletely negligible in the Unitary Gauge we work in. Therefore the Higgsless model results 
for all processes coincide with those in the AIh — > oo limit. 

The cuts in Tab. 1 have been applied either at generation level or as a preliminary step 
to any further analysis. They require containment within the active region of the detectors 
and minimum transverse momentum for all observed partons; a minimum mass separation 
is imposed for all same-family opposite-sign charged leptons, while a minimum separation 
both in mass and in AR is required for jet pairs. At large transverse momentum, jet pairs 
with mass comparable to the mass of electroweak bosons or even larger can merge into 
one single jet when an angular measure like AR{jj) is adopted for reconstructing jets. 
Therefore we have imposed that all partons satisfy AR{jj) > 0.3, a value smaller than 
usually employed in LHC analyses. However, in Sect. 8 we discuss in more detail the effect 
of removing the angular separation constraint or, on the contrary, of imposing a more 
stringent requirement AR{jj) > 0.5. Furthermore, the most forward and most backward 
jets are required to be separated by at least four units in rapidity and their combined 
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mass is forced to be outside the electroweak vector boson mass window; on the contrary 
the mass of the two remaining jets, which we will call central in the following, is required 
to be compatible with the mass of the weak bosons; no jjjijiv) triplet is allowed in the 
neighborhood of the top mass for the £'^£~ + 4j{3£u + 2j) channel. 

Acceptance cuts 

Pt(^^) > 20 GeV ~ 
~ |r?(^^)|<3.0 

PtU) > 30 GeV ~ 
< 6.5 

M{jj) > 60 GeV 

M {£+£-) > 20 GeV 
~M{jfjb) < 70 GeV;M{jfjb) > 100 GeV 

|A??(jfjb)| >4.0 

\M{jjj) - Mtopl > 15 GeV {£+£- + Aj) 
\M{j£u) - Mtopl > 15 GeV {S£u + 2j) 
70 GeV < M{jcjc) < 100 GeV 
ARjjj) > 0.3 

Table 1: Standard acceptance cuts applied in the event generation and present in all results. Here 
j = d, u, s, c, b, g. The last two cuts refer to the £'^£~ + 4j and iv + 4j cases only. 



5. The + 4j channel 

The total cross section for the £'^£~ + Aj channel with the generation cuts in Tab. 1 is 
presented in Tab. 2 as a function of the minimum invariant mass Mcut for the jcjc^'^£~ 
system. These results refer to the mass window 70 GeV < M{jcjc) < 100 GeV and include 
all three perturbative orders. In parentheses the results for the sum of the 0{a%j^j) and 
^{'^EM^D processes. Tab. 2 shows that the cross section is dominated by the 0{a\f^ja^) 
contribution. If we assume a luminosity of L = 200 fb~^ and sum over the electron and 
muon channels the difference between the number of events expected for an infinite mass 
Higgs and a light one is smaller than the expected statistical uncertainty for the 0{a%j^ag) 
processes and no meaningful separation between the two cases can be obtained. 

Therefore, on the generated samples we have applied some additional selection cuts. 
They are shown in Tab. 3 for the i^£~ + 4j channel in the order in which they have 
been implemented. The corresponding distributions are presented in Figs. 7-9 for the SM 
and Higgsless cases. In each figure all previous cuts are applied. The vertical dotted line 
indicates the value of the cut and the arrow indicates which part of the events is kept. For 
the distribution of the largest absolute pseudorapidity of the charged leptons in Fig. 7 and 
the distribution of the smallest transverse momentum of the two central jets in Fig. 8 we 
also present the plots normalized to unit area in order to show more clearly the differences 
in shape between the different cases. 
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Mu = ZUU LreV 


a(fb) 


a(fb) 


a(fb) 


600 


45.8(1.2) 


45.7(1.06) 


45.6(1.01) 


800 


23.3(0.605) 


23.2(0.515) 


23.2(0.471) 


1000 


12.7(0.318) 


12.6(0.263) 


12.6(0.24) 


1200 


7.15(0.177) 


7.1(0.147) 


7.1(0.134) 



Table 2: Total cross section for the + 4j channel after acceptance cuts, Tab. 1, in the mass 

window 70 GeV < M{jcjc) < 100 GeV. In parentheses the resuhs for the 0(a|j,J + 0(a^„a|) 
samples. 
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Figure 7: Top row: distribution of the largest absolute pseudorapidity of the charged leptons 
and the corresponding distribution normalized to unit area. Bottom row;distributions of the mass 
and separation in pseudorapidity of the two tag jets. The numbers refer to the fifi + Aj channel 
only. The set of cuts in Tab. 1 are always imposed. The cuts in Tab. 3 are applied incrementally. 
For instance the distribution of \/S.rj{j f at the bottom of the figure includes the additional cuts 
|?7(^^)| < 2.0 and M{jfji,) > 1000 GeV. Interferences between the different perturbative orders are 
neglected. 
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The essence of the set of cuts in Tab. 3 can be easily understood: they amount to 
requiring that the two tag jets are highly energetic and well separated both from each 
other and from the candidate vector bosons; moreover the two reconstructed vector bosons 
are required to be central and have large transverse momentum. 

In Tab. 4 we present the total cross section with the full set of cuts in Tab. 1 and 
Tab. 3 as the minimum invariant mass for the jcjc^~^^~ system is increased between 600 
and 1200 GeV. These results refer to the mass window 70 GeV < M{jcjc) < 100 GeV. In 
parentheses the results for the 0{a%j^) + 0{a%j^^a1) samples. The last column gives the 
cross sections for the 0{a%i^a1) processes alone; the reported values are computed with 
the Higgs mass taken to infinity, they agree within statistical errors with those obtained 
for Mh = 200 GeV. 

Selection cuts 

|r?(^^)|<2.0 
Mjjfjb) > 1000 GeV 

\^v{jfjb)\ > 4.8 

PtUc) > 60 GeV 

\Av{Vj)\ > 1.1 

AR{£+i~) < 1.0 
PtH^I-) > 200 GeV 
Pi-gjc) > 200 GeV 

Max|r/(j)| > 2.8 

Table 3: Additional cuts for the + 4j channel. 



For both the Higgsless and SILH cases we have computed the probability of what we 
call, for brevity and with a slight abuse of language, the BSM hypothesis at 95% SM exclu- 
sion limit, namely the probability that, assuming the absence of a Higgs boson or that the 
SILH model correctly describes nature, the result of an experimental outcome, with a given 
luminosity, has a chance of less than 5% in the Standard Model (PBSM@95%CL). For this 
we proceed as described in Ref. [30], which we briefly summarize for convenience. We define 
the background B as the expected yield of the 0{a%j^^ag) Z + 4j processes and the signal 
S as the expected number of events from all 0{a%j^) and 0(a^J^^a|) processes, both re- 
stricted to the mass interval 70 GeV < M{jf.jc) < 100 GeV and satisfying all cuts in Tab. 1 
and Tab. 3. These contributions correspond respectively to the representative diagrams of 
Fig. 6 (B) and Figs. 1-5 (S). In other words S corresponds to the events which produce a 
resonant boson peak above the flat background due to V + 4j in the mass distribution of 
the two central jets. B and S are considered as random variables representing the number 
of background and signal events for a possible experimental outcome. B and S are the 
corresponding average values which will be taken equal to the predictions of our simulation 
for a luminosity L = 200 fb~^ and summing over the electron and muon channels. We 
take into account the statistical uncertainty of both B and S assuming a standard Poisson 
distribution with average B and S respectively. The predicted signal cross section is also 
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Figure 8: Top row: distribution of the smallest transverse momentum of the two central jets and 
the corresponding distribution normalized to unit area. Bottom row: distributions of the smallest 
separation in pseudorapidity between the vector bosons and any tag jet and |Ai?| separation between 
charged leptons. The numbers refer to the ^/z + 4j channel only. The set of cuts in Tab. 1 are 
always imposed. The cuts in Tab. 3 are applied incrementally. For instance the distribution of 
|Ai7(yj)| includes all the additional cuts presented in Fig. 7 and in the top row of the present 
figure. Interferences between the two perturbative orders are neglected. 



affected by theoretical uncertainties, so the parameter 5 is itself subject to fluctuations. 
For S we assume, in addition to the statistical fluctuations, a theoretical error defined as 
a flat distribution in the window S ib 30% which, in our opinion, is a reasonable choice to 
account for both pdf's and scale uncertainties for the signal. The processes we are inter- 
ested in require center of mass energies of the order of the TeV and therefore involve rather 
large-x quarks, x « 10^^ 10~^ at a typical scale Q of about 100 GeV. In this region 
the uncertainty due to the parton distribution functions is of the order of 5% [40, 41]. As 
already stated, QCD corrections are in the range of 10% and, as a consequence theoretical 
uncertainties are expected to be well within this order of magnitude. Only statistical fluc- 
tuations have been taken into consideration in the case of B. This is motivated by the fact 
that the background is likely to be well measured experimentally from the region outside 
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Figure 9: Distributions of tlie combined transverse momentum of tlie two central jets and the 
largest absolute pseudorapidity of all jets. The numbers refer to the fifj, + 4j channel only. The set 
of cuts in Tab. 1 are always imposed. The cuts in Tab. 3 are applied incrementally. Interferences 
between the two perturbative orders are neglected. 

the signal peak, so that the theoretical error on Z + 4j is not expected to be an issue at 
the time when real data analysis will be performed. We define the test statistics [39] D 
using the following prescription, 

D = B + S-B (5.1) 

D represents the actual number of events found in the peak taking into account sta- 
tistical fluctuations. The probability distribution of D for the three scenarios is reported 
in Fig. 10 for Mcut = 600 GeV with the full set of cuts in Tab. 1 and Tab. 3 and under 
the constraint 70 GeV < M(jcjc) < 100 GeV. The red curve refers to a Higgs of 200 GeV 
while the green one refers to the SILH model and the blue one to the no-Higgs case. The 
dotted vertical line in the plot marks the 95% exclusion limit for the SM predictions. We 
assume a luminosity L = 200 fb^^ and sum over the and ee final states. 



Mcut 

(GeV) 


no Higg 


s 


SILH 


Mh = 200 GeV 


0(a|,,,a|) 


a(fb) 


PBSM 


a(fb) 


PBSM 


a(fb) 


a(fb) 


600 


0.705(0.154) 


77.1% 


0.625(0.076) 


16.8% 


0.595(0.0484) 


0.022 


800 


0.415(0.102) 


75.3% 


0.358(0.0457) 


16.1% 


0.337(0.025) 


0.010 


1000 


0.209(0.06) 


65.8% 


0.174(0.0245) 


12.3% 


0.162(0.0129) 


0.005 


1200 


0.117(0.0328) 


44.7% 


0.098(0.0141) 


9.91% 


0.092(0.0077) 


0.003 



Table 4: Total cross section for the /i"*"/! + 4j channel, with the full set of cuts in Tab. 1 and 
Tab. 3, as a function of the minimum invariant mass Mcut for the jcjclJ'^ IJ^~ system, in the mass 
window 70 GeV < M{jcjc) < 100 GeV. In parentheses the results for the 0(aL/) + 0{ai^,al) 
samples. The last column gives the cross sections for the 0{a'^,^,ja1) processes alone; the reported 
values are computed with the Higgs mass taken to infinity, they agree within statistical errors with 
those obtained for Mh = 200 GeV. The PBSM probabilities refer to a luminosity of L = 200 fb"^ 
and to the sum of the electron and muon channels. 

As reported in Tab. 4, the probability of an experiment to find a result incompatible 
with the SM at 95%CL assuming that the Higgsless model is realized in Nature is of the 



-13- 



prob. distribution function of 4i[^l cliannel 

LiJO.03 
Q 
Q. 

0.025 

0.02 
0.015 

0.01 
0.005 



-20 20 40 60 80 100 120 

D 

Figure 10: Distributions of the discriminant Eq.(5.1) for the i^£^ + Aj{£ — /i, e) channel for 
L = 200 fb"^ and M^ut = 600 GeV with the fuU set of cuts in Tab. 1 and Tab. 3. The red curve 
refers to a Higgs of 200 GeV while the green one refers to the SILH model and the blue one to 
the no-Higgs case. The dotted vertical line in the plot marks the 95% exclusion limit for the SM 
predictions. 

order of 77% for M^ut = 600 GeV and decreases to about 45% for Mcut = 1200 GeV. At 
first sight this is surprising since the difference between the SM and the Higgsless case 
increases with Mcut if only the ©(al^^) + ©(al^^al) processes are considered. In fact if 
for instance one considers the 0(a|^^) + 0(a|^^a|) ratio of the Higgsless to the SM cross 
sections, this can be seen to be larger at higher Mcut- This behaviour can however be 
understood qualitatively in the following way. We can neglect as a first approximation 
the effect of theory errors which only affects the 0{a%j^j) + C'(a|J^^a|) component. Then, 
the separation, in terms of events for a, L = 200 fb^^ luminosity and two lepton channels, 
between the no~Higgs and the light Higgs case amounts to about 2.5 times the standard 
deviation of the two distributions for Mcut = 600 GeV. For Mcut = 1200 GeV however the 
same separation is about equal to the corresponding standard deviation and therefore the 
overlap between the two distributions is larger and the statistical discriminating power is 
reduced. In other terms, the increase of the relative statistical error as the total number of 
events diminishes at larger Mcut is not compensated by the increase of the ratio between 
the number of genuine scattering events in the Higgsless case compared to the SM one. 

For the SILH model the PBSM@95%CL is only about 17% at most, for Mcut = 
600 GeV. 

For L = 200 fb^^ and summing over the mi and ee final states, the expected total 
rates are about 280/250/240 events for the Higgsless/SILH/SM case for Mcut = 600 GeV 
and 47/40/37 for Mcut = 1200 GeV. The ©(alj^.al) background, in the mass window 
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70 GeV < M{jcjc) < 100 GeV we are focusing on, yields about 220 events for Mcut = 
600 GeV and 33 events for M^ut = 1200 GeV. 

6. The 3£u + 2j channel 



The total cross section for the Siv + 2j channel with the acceptance cuts in Tab. 1 is 
presented in Tab. 5 as a function of the minimum invariant mass Mcut for the 3£u system. 
In parentheses the results for the 0{a%^) processes. 



Mcut 

(GeV) 


no Higgs 


SILH 


Mh = 200 GeV 


a(fb) 


a(fb) 


a(fb) 


600 


0.499(0.242) 


0.462(0.213) 


0.452(0.204) 


800 


0.252(0.134) 


0.227(0.114) 


0.219(0.106) 


1000 


0.139(0.078) 


0.122(0.0635) 


0.117(0.0585) 


1200 


0.0815(0.0476) 


0.071(0.0374) 


0.0675(0.0339) 



Table 5: Total cross section for the 3/ii^ + 2j and 2/zei/ + 2j channels after acceptance cuts, Tab. 1. 
In parentheses the results for the 0{a%^) sample. 



Tab. 5 shows that the cross section for the 0{a%j^ja^) processes is about equal to the 
cross section for the 0{a%j^) contribution. In the following we will consider the full sample 
as our signal since essentially all events contain a Z and a W boson. It is possible to 
improve the discriminating power of the analysis increasing the fraction of 0{a%j^) events 
in the event sample since only those are sensitive to the mechanism of EWSB. Therefore, on 
the generated samples we have applied some additional selection cuts. They are shown in 
Tab. 6 in the order in which they have been implemented. The corresponding distributions 
are presented in Figs. 11, 12. The vertical dotted line indicates the value of the cut and the 
arrow indicates which part of the events is kept. The cuts in Tab. 6, as in the +4:j case, 
force the two tag jets to be highly energetic and well separated and the two (reconstructed) 
vector bosons to be central and to have large transverse momentum. 

Selection cuts 

l^jjfib) > 1000 GeV 
PtH^I-) > 200 GeV 
PTiiJ^) > 200 GeV 
\^r]{jfjb)\ >4.8 
\v{jf,jb)\ > 1-2 
\Ar]iVj)\ > 1.5 
\rj{i^)\ < 2.0 
M{Vj) > 300 GeV 

Table 6: Additional cuts for the Siv + 2j channel. 

The total cross section in attobarns for the Sii' + 2j channel, with the full set of cuts in 
Tab. 1 and Tab. 6, as a function of the minimum invariant mass Mcut is shown in Tab. 7. 
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Figure 11: Top row: distributions of the mass of the two tag jets, the transverse momentum of 
the reconstructed W. Bottom row: distributions of the separation in pseudorapidity of the two tag 
jets and of the largest absolute rapidity of the tag jets. The numbers refer to the /i/x/i + 2j and 
/i/ie + 2j channels only. The set of cuts in Tab. 1 are always imposed. The cuts in Tab. 6 are applied 
incrementally. For instance the \Ari{jfji,)\ at the bottom of the figure includes the additional cuts 
MUfjb) > 1000 GeV andpT(^+^") > 200 GeV, priM > 200 GeV. Interferences between the two 
perturbative orders are neglected. 



In parentheses the results for the 0(a|^^) contribution are reported. The cross section is 
dominated by the 0{a%j^j) contribution. The processes including two QCD vertexes give 
only a small contribution which decreases sharply at larger Mcut- 

In Tab. 7 we also give the PBSM@95%CL for the two BSM scenarios. In the present 
case, since no 0(a|J^^a|) background is present, we use as discriminant S, the sum of the 
events for all 0(a|^^) and 0{a%jy,ja1) processes for a luminosity of L = 200 fb~^, summing 
over the /^/U/x, /i/xe, eee and eefj, channels. We take into account the statistical uncertainty 
assuming a standard Poisson distribution and we assume a theoretical error defined as a 
flat distribution in the window S it 30%. 

The corresponding normalized frequency for the three scenarios is reported in Fig. 13 
for Mcut = 600 GeV. The red curve refers to a Higgs of 200 GeV while the green one refers 
to the SILH model and the blue one to the no-Higgs case. The dotted vertical line in the 
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M(jV) [GeV] 

Figure 12: Distributions of the smallest separation in pseudorapidity between the vector bosons 
and any tag jet, of the smallest absolute pseudorapidity of any charged lepton and of the smallest 
invariant mass of the vector bosons and any tag jet. The numbers refer to the /i/i/i + 2j and 
/i/ie + 2j channels only. The set of cuts in Tab. 1 are always imposed. The cuts in Tab. 6 are 
applied incrementally. Interferences between the two perturbative orders are neglected. 



Mcut 

(GeV) 


no Higg 


s 


SILH 


Mh = 200 GeV 


cr(ab) 


PBSM 


a{ah) 


PBSM 


(j(ab) 


600 


26.9(24.8) 


80.0% 


13.5(11.9) 


19.5% 


8.50(6.95) 


800 


18.8(17.8) 


84.1% 


8.20(7.45) 


24.0% 


4.46(3.72) 


1000 


12.8(12.6) 


72.9% 


5.15(4.79) 


15.9% 


2.21(1.83) 


1200 


8.65(8.55) 


65.7% 


3.08(3.00) 


13.1% 


1.19(1.12) 



Table 7: Total cross section for the 3^iy + 2j and 2^ev + 2j channels in attobarns, with the full set 
of cuts in Tab. 1 and Tab. 6, as a function of the minimum invariant mass Mcut for the iiv system. 
In parentheses the results for the ©(a^^^) sample. The PBSM probabilities refer to a luminosity of 
L ^ 200 fb"^ and to the sum of all electron and muon channels. 



plot marks the 95% exclusion limit for the SM predictions. 

The probability of an experiment to find a result incompatible with the SM at 95%CL 
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prob. distribution function of 2j3lv channel 
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Figure 13: Distributions of the discriminant 13 = 5 for the 3£iy + 2j{£ — ii, e) channel for 
L = 200 fb"^ and M^ut = 600 GeV. The red curve refers to a Higgs of 200 GeV while the green 
one refers to the SILH model and the blue one to the no-Higgs case. The dotted vertical line in 
the plot marks the 95% exclusion limit for the SM predictions. 



assuming that the Higgsless model is realized in Nature is of the order of 80% and decreases 
to about 65% for Mcut = 1200 GeV. Because of the absence of large backgrounds this 
channel has a discriminating power which is in fact higher than the corresponding one for 
£^i~ + 4j final states. However the expected rates are quite small. Only about ten events 
are expected for all combinations of muons and electrons for a luminosity of L = 200 fb~^ 
and Mcut = 600 GeV. It might be extremely difficult to obtain an actual measurement 
for such a tiny predicted rate once all the experimental efficiencies are folded in and all 
sources of isolated leptons from the underlying event or from jets faking isolated leptons 
are taken into account. Clearly this channel would greatly benefit from a larger luminosity. 
The corresponding probabilities for the SILH model vary between 13% and 25%. In the 
3ii^+2j channel the PBSM@95%CL is somewhat larger at Mcut = 800 GeV than at smaller 
value because of the steep decrease of the ©(al^^al) background at larger Mcut- 
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7. The + Aj channel 

In this section we recall the results presented in Ref. [30] and complete them with those 
relative to the SILH model. In addition to the acceptance cuts in Tab. 1 we apply the 
selection cuts in Tab. 8. The corresponding cross sections are shown in Tab. 9. These 
results refer to the mass window 70 GeV < M{jcjc) < 100 GeV and include all three 
perturbative orders. In parentheses the results for the sum of the 0{a%j^^) and 0{a%j^^a1) 
processes which we take as our signal S as for the £'^£~ +4j channel. The last columns gives 
the cross sections for the 0{a'^j^a1) processes alone; the reported values are computed with 
the Higgs mass taken to infinity, they agree within statistical errors with those obtained 
for Mh = 200 GeV. The PBSM probabilities are obtained using the procedure detailed in 
Sect. 5. 

Selection cuts 

\i]{£^)\ < 2.0 
Mjjfjb) > 1000 GeV 
|A7?(j/jb)| >4.8 
PtUc) > 70 GeV 
\ArjiVj)\ > 0.6 
priM > 200 GeV 
PTmiss > 100 GeV 

Table 8: Additional cuts for the iv + Aj channel. 



Mcut 

(GeV) 


no Higg 


s 


SILH 


Mh = 200 GeV 


o(«iM«l) 


a(fb) 


PBSM 


CT(fb) 


PBSM 


a(fb) 


a(fb) 


600 


6.07(1.18) 


96.5% 


5.59(0.704) 


35.9% 


5.41(0.524) 


0.23 


800 


3.76(0.779) 


96.8% 


3.40(0.418) 


29.2% 


3.29(0.309) 


0.13 


1000 


2.26(0.483) 


95.4% 


2.01(0.227) 


19.8% 


1.94(0.169) 


0.08 


1200 


1.32(0.263) 


83.9% 


1.19(0.132) 


16.9% 


1.15(0.094) 


0.05 



Table 9: Total cross section for the fiv + Aj channel, with the full set of cuts in Tab. 1 and Tab. 8, 

as a function of the minimum invariant mass Mcut for the jcjcj^^ system, in tlie mass window 
70 GeV < M{jcjc) < 100 GeV. In parentheses the resuhs for the ©(oL/) + 0{al^al) samples. 
The last columns gives the cross sections for the 0{a%^a1) processes alone; the reported values 
are computed with the Higgs mass taken to infinity, they agree within statistical errors with those 
obtained for Mh = 200 GeV. The PBSM probabilities refer to a luminosity of i = 200 fb"^ and 
to the sum of the electron and muon channels. 

For Mcut = 600 GeV the expected 0{a%j^ja\) background in the mass window 70 GeV < 
M{jcjc) < 100 GeV is about 1950 events with a luminosity of L = 200 fb"''^ and summing 
over the electron and muon channels. Correspondingly, 209 signal events arc expected in 
the light Higgs SM scenario and 474 in the no Higgs case. The corresponding prediction 
in the SILH model is 282 events. 
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Figure 14: Distributions of the discriminant Eq.(5.1) for the £v + 4j(^ — /i, e) channel for L — 
200 fb^^ and M^ut = 600 GeV. The red curve refers to a Higgs of 200 GeV while the green one 
refers to the SILH model and the blue one to the no-Higgs case. The dotted vertical line in the 
plot marks the 95% exclusion limit for the SM predictions. 

The normalized frequency of the discriminant D for the three scenarios is reported in 
Fig. 14 for Mcut = 600 GeV. The red curve refers to a Higgs of 200 GeV while the green 
one refers to the SILH model and the blue one to the no-Higgs case. The dotted vertical 
line in the plot marks the 95% exclusion limit for the SM predictions. 

The iv+Aj channel, as expected, is the one with the best discriminating power between 
the BSM models and the SM. The expected number of events in this final state is about 
ten times larger than in the + 4j channel. 

The probability of an experiment to find a result incompatible with the SM at 95%CL 
assuming that the Higgsless model is realized in Nature is of the order of 96% for Mcut = 
600 GeV and decreases to about 84% for Mcut = 1200 GeV. For the SILH model the 
PBSM@95%CL varies between 35% and 17%. 

8. Combining all channels 

In this section we derive the probability that, assuming that either the Higgsless scenario 
or the instance of SILH model we have considered is realized in Nature, the results of the 
measurements of the + 4j, Mv + 2j and + 4j channels at the LHC yield results 
which are outside the 95% probability region for the SM. In order to do so, it is convenient 
to rephrase the method employed in Ref. [30] and in the previous sections for the case of a 
single measurement. In the new language the generalization to a set of several simultaneous 
measurements will be obvious. Given two models A and B and the probability distributions 
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of some physical observable D, P{D\A) and P{D\B), the 95%CL region for model A can 
be defined considering the probability ratio 



R = 



P{D\A) 



< 1. 



P{D\A) + P{D\B) 

One can consider all the possible D values with R > a and fix a with the condition 



that 



/ 



dD P{D\A) 0{R -a) = 95% 



(8.2) 



For distributions like those in Fig. 10, 13 and 14 this procedure would lead for the 
SM to the one-sided exclusion regions shown by a vertical dotted line. The probability for 
model B to yield a result outside this 95%CL region for A would then be 



PBSMmbfoCL = J dD P{D\B) e{a - R) 



(8.3) 



Eqs.(8.1-8.3) can be easily generalized to a set of observables Di, ....,Dn introducing 
the joint probability P{Di, ....,Dn\A) and transforming the integrals Eqs. (8.2-8.3) to N- 
dimensional ones. In the present case we have assumed the probabilities to be independent 
and defined 



N 



P{Di,....,DN\A) = l[P{Di\A) 



(8.4) 



This method is obviously equivalent to the one, based on the Neyman-Pearson lemma, 
of using the likelihood ratio P{D\A)/P{D\B) for discriminating between two different 

hypotheses for a multidimensional test statistics [39]. 

In Tab. 10 we present the PBSM@95%CL for the Higgsless and SILH model compared 
to the SM predictions for each channel and for their combination for M^ut = 600 GeV. We 
assume L = 200 fb"^ and sum over all electron and muon channels. 



channel 


NOH 


SILH 


3£i/ + 2j 


80.0% 


19.5% 


£+£- +4j 


77.1% 


16.8% 


iu + 4j 


96.5% 


35.9% 


semi-leptonic 


99.4% 


41.2% 


all 


99.9% 


51.5% 



Table 10: PBSM@95%CL for for the individual channels and for their combination with the full set 
of cuts for each channel. The results labeled semi-leptonic refer to the combination of the £~^£~ + Aj 
and Iv + 4j channels. 



Tab. 10 shows, within the limits of the parton level, lowest-order analysis presented in 
this paper, that if no Higgs is present we arc essentially certain that the LHC will obtain 
combined results in the iiy + 4j, + 4j and Siv + 2j channels which will lay outside 
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the 95% CL region of the Standard Model. If however some version of the SILH models 
discussed in Ref. [19] is instead realized the probability of such an outcome drops to about 
50%. 



channel 


NOH 


SILH 




AR = 0.0 


AR = 0.3 


AR = 0.5 


AR = 0.0 


AR = 0.3 


AR = 0.5 


3ii^ + 2j 


80.0% 


80.0% 


80.0% 


19.5% 


19.5% 


19.5% 


£+£- + 4j 


90.1% 


77.1% 


38.5% 


22.5% 


16.8% 


8.5% 


iu + 4j 


99.9% 


96.5% 


79.9% 


51.8% 


35.9% 


25.1% 


semi-leptonic 


99.99% 


99.4% 


90.1% 


58.9% 


41.2% 


26.8% 


all 


99.999% 


99.9% 


98.7% 


66.4% 


51.5% 


39.8% 



Table 11: PBSM@95%CL for the individual channels and for their combination for different values 
of the jet cone separation AR. Both the Higgsless case and the SILH model are considered. 



In Tab. 11 we show how our results depend on the size of the AR separation cut among 
jets, comparing our standard choice AR = 0.3 with the the requirement of a stronger 
separation, AR = 0.5 and with case in which no cone separation among jet is required, 
AR = 0.0. In all cases each jet pair is required to have an invariant mass of at least 
60 GeV. As already mentioned in Sect. 4, at large transverse momentum, jet pairs with 
mass comparable to the mass of electroweak bosons can merge into one single jet when 
an angular measure like AR{jj) is adopted for reconstrTicting jets. Since we insist on 
requiring four jets in the final state these events arc discarded and the cross sections are 
smaller at larger AR. However since the W^s transverse momentum distribution is harder 
in the SILH model and in the Higgsless case than in the SM the PBSM@95%CL decreases 
significantly at larger values of AR for the semileptonic channels. The overall combination 
is less sensitive since the Siu + 2j channel is not affected by the AR cut. 

In several analyses [42] it has been shown that significantly better results in the iden- 
tification of hadronic decays of vector boson, Higgs bosons or supersymmetric partners of 
ordinary particles can be obtained analyzing the substructure of high-pr jets in order to 
separate ordinary QCD jets from those generated by the decay of heavy objects. Such an 
approach might overcome the loss in discrimination between the SM and the BSM scenarios 
at large AR. We leave this possibility for further work. 

9. Conclusions 

We have examined at parton level the processes pp — > i~^£~ + 4j, pp — > Siv + 2j and 
pp ^ iu + Aj including all irreducible backgrounds contributing to these six parton final 
states. We have considered three scenarios: a light Higgs SM framework with Mh = 
200 GeV, one instance of the SILH models and an infinite mass Higgs scenario in order 
to determine whether the two BSM models can be distinguished from the SM at the LHC 
using boson-boson scattering. For the semileptonic channels, the largest background is 
V + 4j at 0{a^j^ag) which can be subtracted looking at the distribution of the invariant 
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mass of the two most central jets in the region outside the weak boson mass window. We 
have estimated the probabihty, in the two BSM scenarios, of finding, combining all set of 
measurements, a result outside the 95% probability range in the Standard Model. This 
probability turns out to be about 99.9% for the Higgsless case and 51.5% for the SILH 
model. These probabilities correspond to an integrated luminosity of L = 200 fb~^ and to 
the sum of all electron and muon channels for a mass of the reconstructed pair of vector 
bosons larger than 600 GcV. Jet resolution plays a crucial role in the present analysis as in 
all processes in which high transverse momentum vector bosons or top particles are present 
and decay hadronically. 
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